Army  Research  Laboratory 


Three-Dimensional  Correlation 
and  Spectral  Functions  for 
Turbulent  Velocities  in 
Homogeneous  and  Surface- 
Blocked  Boundary  Layers 


by  David  Keith  Wilson 


ARL-TR-1 287 


July  1997 


19970728  087 


Approved  for  public  release;  distribution  unlimited. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of 
the  Army  position  unless  so  designated  by  other  authorized  documents. 

Citation  of  manufacturer's  or  trade  names  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  thereof. 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


Army  Research  Laboratory 

Adelphi,MD  20783-1 197 

ARL.TR-1287  ~  July  1997 


Three-Dimensional  Correlation 
and  Spectral  Functions  for 
Turbulent  Velocities  in 
Homogeneous  and  Surface- 
Blocked  Boundary  Layers 


David  Keith  Wilson 

Information  Science  and  Technology  Directorate 


Approved  for  public  release;  distribution  unlimited. 


A  pair  of  three-dimensional  (3D)  models  for  correlation  functions  and  spec¬ 
tra  of  velocity  fluctuations  in  turbulent  boundary  layers  is  presented.  First, 
the  case  of  homogeneous  turbulence  is  considered.  Von  Kar man’s  energy 
spectrum  is  used  to  develop  a  complete  set  of  3D  correlation  and  spectral 
equations.  Second,  it  is  shown  how  the  homogeneous  spectra  can  be  modified 
to  include  the  effect  of  eddy-blocking  at  the  ground.  Assuming  that  the  dis¬ 
turbance  to  the  turbulent  flow  resulting  from  the  blocking  is  irrotational,  an 
equation  is  developed  that  allows  one  to  write  the  vertically  inhomogeneous, 
2D  cross  spectra  as  a  function  of  the  2D  cross  spectra  for  a  homogeneous 
flow.  Although  there  are  only  two  adjustable  parameters  in  the  inhomoge¬ 
neous  model,  the  variance  and  a  length  scale,  the  model  is  shown  to  agree 
quite  well  with  a  variety  of  previous  results  for  the  atmospheric  convective 
boundary  layer  (CBL). 
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1.  Introduction 


Three-dimensional  (3D)  spectral  and  correlation  models  of  turbulence  in  at¬ 
mospheric  boundary  layers  are  required  in  a  variety  of  applications,  such 
as  electromagnetic  and  acoustic  wave  propagation,  turbulent  transport  and 
dispersion,  response  of  structures  exposed  to  the  wind,  and  large-eddy  sim¬ 
ulation.  In  particular,  the  work  described  in  this  report  is  motivated  by  the 
need  for  good  spectral  models  of  large-scale  turbulence  in  acoustical  propa¬ 
gation  applications. 

It  is  not  difficult  to  develop  good  3D  spectral  and  correlation  models  if 
one  is  interested  only  in  small-scale  inertial  subrange  turbulence  (fig.  1). 
Owing  largely  to  the  work  of  Kolmogorov  (1941),  the  spectral  properties 
of  the  inertial  subrange  are  well  understood.  The  eddies  belonging  to  the 
inertial  subrange  have  spatial  scales  smaller  than  the  instability  generating 
the  turbulence,  and  larger  than  the  scales  at  which  the  turbulent  kinetic 
energy  is  dissipated  by  viscosity.  The  large  eddies  generated  directly  by  the 
instability  belong  to  the  energy-containing  subrange,  and  the  very  small, 
dissipating  eddies  belong  to  the  dissipation  subrange. 

For  an  atmospheric  convective  boundary  layer  (CBL),  the  energy  input  oc¬ 
curs  at  scales  on  the  order  of  Zj,  the  inversion  layer  height,  which  is  typically 
in  the  range  from  0.5  to  2  km.  Dissipation  occurs  on  scales  ~  1  mm.  We  see 
that  there  is  an  extremely  broad  inertial  subrange  in  atmospheric  turbulence, 
spanning  five  to  six  orders  of  spatial  magnitude.  Nonetheless,  if  turbulent 


Figure  1:  Three  subranges  of  high-Reynolds-number  turbulence,  (a)  Turbu¬ 
lent  kinetic  energy  (TKE)  spectrum,  KE{K)/ff^.  (b)  TKE  production  spec¬ 
trum,  KdT{K)/dK.  (App  A  gives  actual  equations  plotted  here.) 
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scales  in  either  the  energy-containing  or  dissipation  subranges  play  a  role  in 
a  given  application,  or  if  interactions  between  the  eddies  and  ground  surface 
(called  blocking  effects)  are  important,  inertial  subrange  spectral  modeling 
will  fail. 

In  outdoor  acoustic  propagation,  the  wavelengths  of  interest  can  be  as  long  as 
30  m,  and  heights  all  the  way  down  to  ground  level  are  important.  Turbulent 
eddies  having  dimensions  on  the  order  of  or  greater  than  the  wavelength  pri¬ 
marily  determine  the  scattered  field.  Hence  classical  inertial  subrange  scaling 
cannot  be  used  in  many  acoustical  applications,  and  models  including  the 
energy-containing  subrange  and  ground  blocking  effects  are  required.  Turbu¬ 
lent  structure  in  the  dissipation  subrange  is  normally  unimportant,  though. 

Despite  the  importance  of  having  good  3D  statistical  models  for  atmospheric 
turbulence,  previous  experimental  and  modeling  efforts  have  largely  focused 
on  single-point  variances  and  ID  spectra  taken  in  the  direction  of  the  mean 
wind.  This  is  because  such  quantities  can  be  determined  from  stationary 
towers.  For  ID  spectra,  Taylor’s  frozen  turbulence  hypothesis  is  used  to 
convert  a  time  series  recorded  at  the  stationary  tower  to  a  spatial  series. 
(Discussions  of  ID  spectral  measurements  and  models  are  found  in  Kaimal 
et  ai  (1972),  (1976),  H0jstrup  (1982),  and  Panofsky  and  Dutton  (1984).) 
Although  single-point  variances  and  ID  spectral  models  provide  some  basis 
and  constraints  on  which  to  build  a  multidimensional  model,  they  are  not 
sufficient  to  fully  define  one. 

In  principle,  by  assuming  homogeneity  and  isotropy,  one  could  actually  derive 
multidimensional  spectra  and  correlation  models  solely  from  ID  turbulence 
measurements.  The  necessary  relationships  have  been  known  for  decades 
(Batchelor,  1953).  Unfortunately,  one  cannot  apply  these  relationships  di¬ 
rectly:  atmospheric  turbulence  is  known  to  be  highly  inhomogeneous  and 
anisotropic,  largely  because  of  interactions  with  the  ground. 

Despite  the  inherent  difficulties,  some  notable  efforts  have  been  made  to  de¬ 
velop  multidimensional  spectral/correlation  models  suitable  for  atmospheric 
turbulence.  Kristensen  et  al  (1989)  developed  a  general  3D  model  for  tur¬ 
bulent  velocity  spectra  in  homogeneous,  anisotropic  boundary  layers.  Mann 
(1994)  considered  3D  spectra  for  turbulence  generated  by  shear  instabilities, 
allowing  vertical  inhomogeneities  in  his  model.  The  work  in  this  report  builds 
on  a  spectral  theory  developed  by  Hunt  and  Graham  (1978),  who  considered 
the  effect  of  a  boundary  on  turbulence  in  the  absence  of  mean  flow.  As  a 
result,  the  model  here  is  intended  for  the  atmospheric  CBL.  By  this  I  mean 
that  turbulence  generation  by  convective  instabilities  (heating  by  the  ground 
surface  and  the  resulting  density  contrasts  in  the  overlying  air)  dominates 
generation  by  shear  instabilities  (the  wind  blowing  over  the  ground). 

One  of  the  principal  benefits  of  building  on  Hunt  and  Graham’s  (1978)  results 
is  that  the  resulting  spectral  equations  give  the  inhomogeneous,  boundary- 
blocked  spectra  in  terms  of  the  spectra  for  homogeneous  turbulence.  Since 
the  homogeneous  case  can  be  analyzed  fully,  as  mentioned  above,  a  complete 
spectral  model  results. 


This  report  is  divided  into  three  main  sections.  Section  2  covers  the  case 
of  homogeneous,  isotropic  turbulence.  I  provide  some  tutorial  background 
for  readers  unfamiliar  with  this  subject,  and  I  develop  a  complete  set  of 
spectral/correlation  equations  for  the  homogeneous  case.  The  main  purpose 
of  the  section,  however,  is  to  derive  the  spectral  results  for  homogeneous 
turbulence  needed  by  the  inhomogeneous,  blocked  model. 

The  inhomogeneous,  blocked  model  is  developed  in  section  3.  I  derive  a 
general  equation  describing  the  blocking  on  2D  cross  spectra,  by  assuming 
that  the  boundary  does  not  modify  the  vorticity  field.  All  other  correlations 
and  spectra  can  be  determined  from  the  2D  cross  spectra  by  calculation  of 
Fourier  transforms. 

Last,  in  section  4,  I  compare  the  inhomogeneous,  blocked  model  to  several 
experimental  and  computationally  generated  data  sets,  as  well  as  to  some 
previously  developed  models.  Of  course,  as  was  discussed  above,  it  is  diffi¬ 
cult  to  perform  multidimensional  correlation/ spectral  measurements  in  the 
atmosphere.  However,  quantities  that  are  more  easily  measured,  such  as  vari¬ 
ances,  integral  length  scales,  and  ID  spectra,  are  special  cases  of  the  full  3D 
correlation/spectral  model,  and  I  expect  the  model  described  here  to  behave 
well  for  these  special,  limiting  cases.  Hence  it  is  quite  worthwhile  to  test  the 
model  with  measurements  of  these  simpler  quantities. 
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2.  Correlation/Spectral  Model  for  Homogeneous, 
Isotropic  Turbulence 


2.1  Preliminaries 

Suppose  we  measure  turbulent  velocity  fluctuations  at  two  points  in  space. 
Let  us  indicate  these  two  points  as  x  =  (x^y^z)  and  x'  =  (x',  The 

coordinate  system  is  oriented  so  that  2:  is  the  vertical  direction.  Orientation 
of  the  horizontal  axes  is  arbitrary,  since  I  am  not  considering  turbulence 
production  by  a  mean  horizontal  wind  shear  in  this  report.  The  correlation 
function  between  the  measured  velocities  is 

Rij{x,x')  =  {ui{x)uj{x')),  (2.1) 

where  the  angle  brackets  indicate  ensemble  averaging,  and  ui  and  uj  are 
the  velocity  fluctuations  at  the  two  points,  oriented  along  the  ith  and  jth 
coordinate  axes  (where  x\  is  the  same  as  x,  X2  =  y,  and  X3  =  2:). 

By  definition,  in  homogeneous  turbulence,  the  correlation  function  depends 
only  on  the  spatial  separation  between  x  and  x'.  Hence  we  may  write 

Rij{r)  =  (tii(x)uj(x')),  (2.2) 


where  r  =  x'  —  x. 

The  spectra  are  Fourier  transforms  of  the  correlation  function.  The  Fourier 
transform  convention  used  in  this  report  is 

/(«)  =  ^  fir)  exp(-mr)  dr  (2.3) 


for  the  forward  transform,  and 


/oo  ^ 

/(K)exp(iKr)  dK. 

-OO 


(2.4) 


for  the  inverse  transform,  in  which  k  is  the  wavenumber.  For  example,  the 
ID  cross  spectrum  is 

1 

Qij{t^i',r2,r3)  =  —  J  Rij{ri,r2,r3)exp{-iKiri)dri.  (2.5) 


By  cross  spectrum,  I  mean  a  spectrum  involving  two  spatially  separated 
points.  The  ordinary,  single-point,  ID  spectrum  would  be  0ij(Ki;0, 0).  A 
cross  spectrum  can  be  thought  of  as  a  mixed  spatial  correlation/spectral 
function. 

It  is  only  a  matter  of  convention  that  the  Fourier  transform  was  taken  in  the 
xi-direction  (i.e.,  with  respect  to  ri)  in  equation  (2.5).  Because  of  isotropy, 
we  could  have  used  either  the  X2-  or  xa-axis  instead,  and  still  had  a  complete 
set  of  ID  cross  spectra.  For  example,  consider  the  ID  cross  spectrum  with 
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Fourier  transform  in  the  X3-direction  ©22(’'x»  ^||  1 Let  us  use  the  customary 
right-handed  coordinate  system,  with  the  xi-axis  coming  out  of  the  page, 
X2  to  the  right,  and  X3  upwards.  Looking  in  the  direction  of  the  velocity 
component  (x2-axis),  we  see  the  displacement  to  the  right,  the  displace¬ 
ment  ry  to  the  front,  and  the  Fourier  transform  axis  upwards.  By  reorienting 
ourselves  so  that  we  are  looking  down  the  X3-axis  with  the  X2-axis  to  the 
right,  we  find  the  relationship  ©33(k;’"x>^|1)  =  ©22(’”X5^’1|)'^)-  Many  other 
relationships  follow  similarly. 

Besides  the  ID  cross  spectra,  I  define  2D  cross  spectra  (f>ij{Ki,  K2-,r3),  and 
3D  spectra  (ki,  «2)  as  follows: 


1 

(ki,  K2;  r3)  =  -^  J  ^2,  ^3)  exp(-iKiri  -  iK2r2)  dn  dr2 

1 

=  TT  I  ©ij(«;i;^2,r3)exp(-i/C2r2)dr-2, 

27r  J —00 


(2.6) 


#y  (Kl,K2,  K3) 


1 

— 5-  /  Rij(ri,r2,r3)exp{-iKiri  —  iK2r2  —  iK3r3)dridr2dr3 

7—00 

^  f  (f>ij{Ki,K2]r3)exp{-iK3r3)dr3.  (2.7) 

iTt  J-00 


The  eventual  goal  in  these  modeling  efforts  is  to  determine  Rij,  ©ijj  <Pijf 
.  Of  course,  because  of  the  Fourier  transform  interrelationships,  only  one 

i'J 

of  these  functions  is  independent. 

There  are  other  useful  symmetry  relationships.  Switching  Ui  and  uj  in  the 
definition  of  the  correlation  function,  and  applying  homogeneity,  we  find 

Rij{ri,  r2,  r3)  =  Rji{-ri,  -r2,  -rs)-  (2-8) 

Making  use  of  this  identity  and  the  fact  that  Rij  is  a  real  function,  we  find, 
from  the  Fourier  transform  definitions, 

©y(Ki,  K2,  r3)  =  ©ji(-Ki;  -r2,  -r-3)  =  -r2,  -rs),  (2.9) 

<?i'i7(«l,K2;r'3)  =  (l>ji{-Kl,-K2]-r3)  =  (f>ji{Kl,K2;-r3), 

(2.10) 


and 


#y(Ki,K2,K3)  =  ^ji(-«l,-K2,-K3)  =  $ji(Kl,«2,«3)- 

(2.11) 


2.2  Energy  and  Three-Dimensional  Spectra 

The  starting  point  in  developing  correlations  and  spectra  for  homogeneous, 
isotropic  turbulence  is  the  3D  turbulent  kinetic  energy  (TKE)  spectral  den¬ 
sity  E(«),  where  k  =  (/ci,  K2,  «3),  and  k  =  |ic|.  (In  the  turbulence  literature. 
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TKE  usually  implicitly  refers  to  the  TKE  per  unit  mass.  I  adhere  to  that 
convention  in  this  report.)  E{k)  is  defined  as  the  TKE  per  unit  wavenumber 
magnitude.  If  one  thinks  of  the  TKE  in  3D  wavenumber  space,  E{k)  is  the 
energy  in  a  shell  of  radius  k.  Since  the  TKE  spectral  density  at  a  given  point 
in  wavenumber  space  is  $ii(K)/2  (summation  over  repeated  roman  indices 
is  implied),  and  the  total  energy  in  a  shell  of  radius  k  is  times  that 
amount, 

E{k)  =  2TlK^^ii{K).  (2.12) 

The  total  kinetic  energy  in  the  flow  is 

roo  2  roo  roo  roo 

E=  E{K)dK=-  /  /  /  ^>u(k)  dKi 

Jo  2i  J —  CO  J —  OO  j —  OO  /'O  1 


The  success  of  the  correlation/spectral  model  depends  on  making  a  good 
choice  for  E{k),  The  main  criteria  are  (1)  it  must  agree  well  with  data:  and 
(2)  it  must  be  analytically  convenient  to  manipulate.  The  following  form  for 
the  3D  energy  spectrum  satisfies  both  criteria  reasonably  well  within  the 
energy-containing  and  inertial  subranges  (see  app  A): 


_  4r{u  +  5/2) 

V^r(z/)  (1  +  k2£2)^+5/2’ 


(2.14) 


where  is  the  variance,  k  the  radial  wavenumber,  £  the  length  scale,  T 
the  gamma  function,  and  u  sets  the  power  law  dependence  in  the  inertial 
subrange  ':$>  1).  With  appropriate  choices  for  <7^,  and  z/,  equation  (2.14) 
is  equivalent  to  von  Karman’s  (1948)  model  for  E{k).  In  particular,  setting 
z/  to  1/3  results  in  Kolmogorov’s  “-5/3  power  law  for  the  inertial  subrange 
{k£  1).  I  will  not  explicitly  set  u  to  1/3  until  we  obtain  final  results,  since 

to  do  so  would  sacrifice  generality  without  appreciably  simplifying  any  of 
the  spectral  expressions. 

Given  the  energy  spectrum,  it  is  a  simple  matter  to  determine  the  3D  spectra 
One  uses  the  following  well-known  formula,  derived,  for  example,  in 
Batchelor  (1953): 


^ij{K) 


E(k) 


KiKj 


)■ 


(2.15) 


Longitudinal  Correlations  and  Spectra 

As  a  next  step,  we  could  integrate  the  3D  spectra  to  find  the  2D  cross 
spectra  (pij.  It  is  somewhat  more  instructive,  however,  to  first  consider  ID 
spectra.  Given  equation  (2.14),  and  the  assumption  of  homogeneous,  isotropic 
turbulence,  it  turns  out  that  ID  spectral  densities  and  correlation  functions 
can  be  derived  directly. 

Let  us  begin  with  the  ID  longitudinal  autospectrum.  By  longitudinal,  I  mean 
that  the  velocity  component  is  parallel  to  the  direction  of  the  wavenumber 


axis:  i.e.,  eii(«:,  0, 0),  622(0,  ac,  0),  or  633(0, 0,  k).  Note  that  isotropy  implies 
that  these  three  spectra  are  the  same.  Let  us  define  a  normalized,  longitu¬ 
dinal  spectrum  as  /{k)  =  6ii(/^, 0, 0)/cr^,  which  can  be  shown  to  be  related 
to  the  energy  spectrum  as  follows  (Batchelor,  1953): 


E{k) 


2  3 
=  a  K 


dK 


1  df{K) 
K  dK 


(2.16) 


By  integrating  equation  (2.16),  it  is  straightforward  to  show  that  the  nor¬ 
malized,  longitudinal  spectrum  is 


f(  ^  -  r(^  +  V2)  e 

V^T{u)  (H-k2£2)./+1/2- 


(2.17) 


The  longitudinal  spectrum  is  plotted  as  the  solid  line  in  figure  2,  for  i/  —  1/3. 

We  can  derive  the  longitudinal  correlation  function  cr^/(r)  by  taking  the  in¬ 
verse  Fourier  transform  of  equation  (2.17).  Using  integral  (3.771.1)  in  Grad- 
shteyn  and  Ryzhik’s  (1994)  tables,  and  the  fact  that  the  autospectrum  is  an 
even  function,  we  find  the  normalized  longitudinal  correlation  to  be 


(2.18) 


where  is  the  modified  Bessel  function  of  the  second  kind  of  order  u. 

Since  modified  Bessel  functions  play  an  important  role  throughout  the  re¬ 
mainder  of  this  report,  it  is  worth  taking  a  moment  to  discuss  a  couple  of 


Figure  2:  One-dimensional  autospectral  density  model  for  i/  =  1/3. 
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their  properties.  For  small  arguments  |  1  and  v  <2, 


Km 


2u  \2j  2{l-u)\2j 


The  leftmost  of  the  three  terms  on  the  left  is  always  the  leading  term  in 
the  series;  the  next  most  significant  term  in  the  series  is  either  the  second 
or  third,  depending  on  the  particular  value  of  u.  Substituting  the  first  term 
into  equation  (2.18).  it  is  eaisy  to  show  that  /(O)  =  1,  as  it  must.  For  large 
arguments  ^  ^  1,  becomes  a  decaying  exponential: 

~  y^exp(-0.  (2.20) 


Equation  (2.18)  is  plotted  as  the  solid  line  in  figure  3.  It  must  be  kept 
clear  that  /(r)  is  the  normalized  autocorrelation  only  when  the  displace¬ 
ment  is  in  the  longitudinal  direction.  That  is,  i?ii(r,  0, 0)  =  Il22(0,  r,  0)  = 
.^33(0, 0,r)  =  cr^/(r).  It  is  not  generally  true,  even  for  an  isotropic  vector 
field,  that  i?ii(r,0,0)  equals  iiii(0,r,0).  The  normalized  transverse  correla¬ 
tion  function  g{r)  is  addressed  in  the  next  subsection. 


2.4  Transverse  Correlations  and  Spectra 

So  far,  we  have  examined  the  autocorrelation  and  autospectral  functions 
when  the  displacement  (or  the  wavenumber  component)  is  parallel  to  the 
direction  of  the  velocity  component.  In  order  to  study  more  general  types  of 
correlations,  we  need  to  develop  an  expression  for  the  complementary  case. 


Figure  3;  Autocorrelation  function  model  for  u  =  1/3. 


where  the  displacement  is  perpendicular  to  the  velocity  component.  This  is 
readily  accomplished  with  the  theory  of  isotropic  vector  fields  in  multiple 
dimensions,  such  as  discussed  in  Batchelor  (1953).  For  incompressible  fiow 
{dui/dxi  =  0),  the  transverse  correlation  g{r)  can  be  computed  from  the 
longitudinal  /(r)  with  the  following  formula  (Batchelor,  1953): 

g(^r)^  f{r)  +  ~.  (2.21) 


In  this  case,  using  equation  (2.18),  we  find 

1  /r\'^ri^  +  2. 


=  ¥kr)(j) 


2^r(i/) 


(2.22) 


(The  argument  of  the  modified  Bessel  functions  is  implicitly  r/£,  unless  oth¬ 
erwise  specified.)  In  deriving  the  second  and  third  forms  above  for  g{r),  I 
made  use  of  the  differentiation  formula  for  modified  Bessel  functions 

KiO  =  -I  [Ku-iiO  +  ^.+i(0] ,  (2-23) 


and  the  recurrence  relation 

=  +  (2-24) 


Most  generally,  an  isotropic,  homogeneous  vector  field  has  a  correlation  ten¬ 
sor  of  the  form  (Batchelor,  1953) 


(2.25) 


where  n  is  the  spatial  separation  along  the  ith  axis,  and  +  r^. 

Using  the  final  form  in  equation  (2.22)  for  the  transverse  correlation,  we 
obtain  the  full  isotropic  correlation  tensor: 


A,V)  =  54m  ^  -  ?^)  [O'  +  5  (9  }  • 


(2.26) 


Equation  (2.26)  is  the  complete  correlation  model  for  the  velocity  fluctua¬ 
tions  in  homogeneous,  isotropic  turbulence.  It  is  plotted  in  figure  3  for  i  =  j 
and  various  angles  between  the  displacement  vector  and  the  velocity.  Note 
that  the  transverse  correlation  {i  =  j,  with  rj  =  rj  =  0)  actually  becomes 
negative  for  large  values  of  r/£.  Although  this  behavior  may  seem  surprising 
initially,  it  is  a  physically  necessary  property  of  the  transverse  correlation 
(Batchelor,  1953),  not  just  an  artifact  of  using  the  von  Karman  energy  spec¬ 
trum. 
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Finally,  consider  the  normalized  transverse  spectral  density  function,  which 
is  the  Fourier  transform  of  g{r).  In  the  final  form  in  equation  (2.22),  there  are 
two  terms,  one  involving  and  the  other  The  Fourier  transform  of 

each  term  follows  from  (6.726.4)  in  Gradshteyn  and  Ryzhik’s  (1994)  tables. 
After  some  algebra,  we  find 

r(i/  +  i/2)  e  (  i/  +  i/2\ 

^Yiv)  +  +  ^  l  +  K^ej-  (2.27) 

To  check  equation  (2.27),  consider  the  inertial  subrange  limit  k\(.  1.  Then 
the  term  in  brackets  becomes  r/  +  1  =  4/3.  Comparison  to  equation  (2.17) 
shows  that  gi^n)  =  (4/3)/(/t)  in  the  inertial  subrange,  a  well-known  result 
from  the  theory  of  isotropic  vector  fields  (Panofsky  and  Dutton,  1984). 


2.5  One-Dimensional  Cross  Spectra 

We  now  have  determined  all  the  ID  autospectra,  in  free  space.  For  exam¬ 
ple,  ©ii(k,0,0)  =  ©22(0,  K,0)  =  and  ©ii(0,  k:,0)  =  ©22(«^5  0,0)  = 

The  other  cases  can  be  found  by  rotation  of  the  coordinate  axes.  The 
next  step  is  to  find  the  ID  cross  spectra,  which  we  determine  by  transform¬ 
ing  the  correlation  function  (2.26).  For  concreteness,  I  take  the  direction  of 
the  Fourier  transform  to  be  the  Xi  direction,  and  indicate  the  result  using 
the  following  notation: 

1 

©ii(^i;r2,r3)  =  —  /  Rij{rur2,rs)exp{iKiri)  dri. 

/tt  J-oo  ^2.28) 


Consider  first  only  the  longitudinal  part  of  the  spectrum,  the  first  term  in 
equation  (2.26).  Defining  +  rf,  we  have  the  integral 


7r2^~^r{i/) 


r( 


2  ^  ^2 


X 


cos(«:iri)  dr\. 


(2.29) 


This  integral  can  be  calculated  in  closed  form;  it  is  given  as  (6.726.4)  in 
Gradshteyn  and  Rhyzhik’s  (1994)  tables.  We  find 


i  f  ^  \i/+l/2 

/(«i;r2,r3)  =  ^2''-i/2r(t/)  Vl  +  K?^V  -^^^<^+1/2  (0  , 


(2.30) 


where  =  (r^  -|-  r3)(l  + 

The  computation  becomes  somewhat  more  complicated  when  we  include  the 
transverse  correlation,  although  results  still  can  be  obtained  in  closed  form. 
Writing  out  the  correlation  function  Rn  as  given  by  equation  (2.26),  and 
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simplifying  with  the  recurrence  formula  (eq  (2.24)),  we  find 


Rn{ri,r2,r3) 


2'^-ir(i/) 


r ,  ( r 


-1 

K^-i 


(2.31) 


The  Fourier  transforms  of  both  terms  in  the  square  brackets  are  obtained 
readily.  The  first  term  integrates  in  the  same  manner  as  the  longitudinal 
spectrum,  whereas  the  second  term  has  the  same  form  if  one  replaces  v  by 
v-l.  The  result  of  the  integration  is 


©u(Ki;7'2,r3)  = 


aH 


0F2‘^-V2r(jy) 


\  t/+l/2  r 


(0  -  iO 


(2.32) 


Equation  (2.31)  is  also  valid  for  i?22,  if  we  set  r\  =  r?  +  r|.  Calculation 
of  the  spectrum  ©22(«i,^2,^3))  however,  is  somewhat  more  difficult  than 
that  of  ©11,  because  now  r±_  depends  on  the  variable  of  integration  in  the 
Fourier  transform,  ri.  It  turns  out  that  we  can  use  the  recurrence  relation  (eq 
(2.24))  to  put  the  correlation  function  in  a  form  that  is  more  conveniently 
integrated: 


i?22(ri,r2,r3)  =  G)  ”  \  G)  ^  (^) 


(2.33) 


Each  of  the  terms  in  square  brackets  now  is  readily  transformed  with  Grad- 
shteyn  and  Rhyzhik’s  equation  (6.726.4),  with  the  result 


©22(«;i;r2,r3) 


V^2^-i/2r(z.)  Vl  +  «?^V 

£  r^(l  +  k?£^) 

(i/  +  1)K^+i/2  (0  -  2(1  +  (0  +  ^  ^^-1/2  (0 


(2.34) 


Because  of  isotropy,  ©33(Ki,r2,r3)  must  also  be  given  by  equation  (2.34), 
except  with  the  indices  2  and  3  interchanged. 

The  correlation  function  R12  follows  from  equation  (2.26)  and  the  recurrence 
relation  (eq  2.24)): 

Bi2(n,r2.r3)  =  j^(0  (2.35) 

The  Fourier  transform  with  respect  to  ri  can  be  determined  with  the  help 
of  Gradshteyn  and  Rhyzhik’s  equation  (6.726.3),  the  result  being 


©i2('«i;^2,»'3) 


ah2Kie  (  i 
0F2‘^-i/2r(i/)  \i  +  Kle^J 


11 


The  cross  spectrum  ©13  is  identical  to  ©12,  except  that  in  the  equation 
above  is  replaced  by  r^. 

The  correlation  function  i?23  is  the  same  as  R\2,  except  with  rs  replacing  ri- 
Taking  the  Fomier  transform  and  using  Gradshteyn  and  Rhyzhik’s  equation 
(6.726.4),  one  finds 


023(Ki;r2,r3)  = 


-1/2 


Kv-ll2  {0  ■ 


(2.37) 


When  dealing  with  spectra,  it  is  important  to  keep  clear  the  distinction  be¬ 
tween  their  one-  and  two-sided  versions  (Bendat  and  Piersol,  1986).  The 
functions  ©^  discussed  above  are  two-sided  spectra,  being  defined  for  both 
positive  and  negative  ki-  In  fact,  Qij{—Ki,0,z,z)  —  Qji{Ki,0,z,z).  Experi¬ 
mentalists  more  commonly  measure  one-sided  spectra  Faa-,  which  are  defined 
as  zero  for  negative  /ti  and  twice  Qaa  for  positive  ki  .  Note  that  the  variance 
is  recovered  by  integration  of  the  two-sided  spectral  density  from  kj  =  —  00 
to  -t-00,  whereas  the  variance  is  recovered  from  the  one-sided  version  by 
integration  from  0  to  -t-00. 


2.6  Two-Dimensional  Cross  Spectra 

We  can  find  the  2D  cross  spectra  by  calculating  either  forward  Fourier  trans¬ 
forms  of  the  ID  cross  spectra,  or  inverse  transforms  of  the  3D  spectra.  It  is 
somewhat  easier  to  use  the  latter  approach.  Solutions  to  the  necessary  inte¬ 
grals  can  be  found  in  Gradshteyn  and  Rhyzhik  (1994),  equations  (3.771.2) 
and  (3.771.5).  We  find 


<pn{Ki,  K2;rz)  = 


V+1 


7r2i/r(i/)(l  -I- 


(2.38) 


<^33(1^1,  «2)  T’s)  —  ^2i'-t-lr(i/)(l  +  (C/i)  5  (2.39) 


^  ^ . \ 

Mki,  K2,r3)  +  ^2py+2^>'+‘2  (C*)  > 


(2.40) 


013(«1)K2;^3)  = 


ia^Kie^Ck+^ 

7r2-+T(l/)(l  +  4£2)t.+3/2^‘'+l 


where  Cft  =  (ra/^)^!  +  xhe  remaining  spectra  follow  from  horizontal 
isotropy: 


and 


<f>2z{K^2^X'l]f'z)  =  <?^13(«1,«2;^3)- 


Frequently  we  are  most  interested  in  the  special  case  of  2D  spectra  having 
no  vertical  separation.  To  evaluate  this  limit,  we  can  use  the  expansion  of 
the  Bessel  function  for  small  equation  (2.19).  Substituting  the  expansion 
into  equations  (2.38)  to  (2.41),  we  find 


^ii(ki,«2;0) 


7r(l  + 


5  +  (»"  +  1) 


IT^J’ 


(2.42) 


^33(ki,«2;0) 


v{y  +  l)(T^fc^£‘* 

7r(l  +  4^2)..+2  ’ 


(2.43) 


^12(ki,K2;0)  =  - 


v{y  +  \)a'^K\K2^^ 

7r(l  +  ’ 


(2.44) 


^I3(«i5  K2;  0)  =  0.  (2.45) 

In  deriving  equations  (2.42)  to  (2.45),  use  was  made  of  the  relationship 

r(i/  +  i)  =  i^r(i/). 
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3.  Inhomogeneous  CBL  Model 


All  classes  of  spectra  and  correlation  functions  for  homogeneous,  isotropic 
turbulence  were  derived  in  the  previous  section.  In  this  section,  I  introduce 
the  complication  of  a  solid,  free  slip  boundary  at  2  =  0.  The  main  effect  is  to 
block  the  large  eddies,  forcing  their  vertical  velocity  to  zero  at  the  boundary. 
This  introduces  vertical  inhomogeneity  into  the  correlations  and  spectra. 


3.1  Modification  of  Two-Dimensional  Cross  Spectra  Near 
a  Boundary 

Let  us  decompose  the  turbulent  velocity  fluctuations  Ui  into  contributions 
from  the  homogeneous  flow  (H),  such  as  would  exist  in  the  absence  of  the 
boundary,  and  a  contribution  from  blocking  at  the  boundary  (B): 

Ui{x,y,z)  =  u\^\x,y,  z)  +  u\^\x,y,  z).  (3.1) 

Since  1x3  must  vanish  at  the  boundary,  this  solution  is  subject  to  the  condition 

4^^  {x,  y,  0)  +  4^^  (^)  2/)  0)  =  0.  (3.2) 

Because  I  am  not  considering  mean  shear,  and  viscous  effects  are  negligible 
for  the  high  Reynolds  numbers  characteristic  of  the  atmosphere,  the  bound¬ 
ary  does  not  modify  the  vorticity  field  (Hunt  and  Graham,  1978).  Hence,  the 
boundary  contribution  to  the  velocity  field  can  be  modeled  by  irrotational, 
inviscid  flow  theory,  and  we  can  define  a  velocity  potential  ^  such  that 

uf'>  =  V$,  (3.3) 

where  $  satisfies  Laplace’s  equation: 

=  0.  (3.4) 

Taking  the  2D  Fourier  transform  of  equation  (3.4),  we  have 

(-4  +  ^  j  «2,  z)  =  0,  (3.5) 

where  =  /cf  +  /C2)  and  the  “hat”  indicates  the  2D  Fourier  transform.  The 
solution  to  this  equation,  subject  to  the  condition  that  the  influence  of  the 
boundary  becomes  negligible  as  z  — >  oo,  is 

$  =  Aexp{—Khz).  (3.6) 

Prom  the  boundary  condition  (eq  (3.2)), 

A  =  uf\Ki,K2,0)/Kh.  (3.7) 
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(3.8) 


By  Fourier  transforming  equation  (3.3),  we  can  now  easily  show  that 

u\^\ki,K2,z)  =  mi{Ki,  K2)  exp{-Khz)ui^\Ki,  K2,0), 


in  which 

mi{Ki,K2)  =  (iKi/Kh,iK,2/x-h,—^)-  (2-^) 

(B) 

We  have  now  determined  the  effect  of  the  boundary  by  solving  for  in 
terms  of  In  principle,  this  allows  us  to  determine  all  the  correlations 
and  spectra.  The  2D  cross  spectra  are  equal  to 

4>ij{Kl,K2,Z,z')  =  {UiiKl,K2,z)u*{Kl,K2,z')).  (3.10) 

Since,  from  equations  (3.1)  and  (3.8), 

Uj(/Ci,  K2,  z)  =  K2,  z)  +  mi  exp(-/c/i2)4^^(Ki,  K2, 0), 


we  have  the  result 

<^y(Kl,  K2;  z')  =  (j)^^\ni,K2]  Z'  -  Z) 

+  e'"'''‘^'mj(/ti,  K2)4f  ^*('^1'  '^2;  z)  ^2) 

+  (ki,  K2)4f  '^2;  z') 

+  K2)mj(Ki,  K2)(^if\Ki,  K2;0). 

The  second  term  on  the  right  was  simplified  by 

<pjp(Ki,K2;-z)  =  <^jP*(Ki,K2;z),  (3.13) 

which  can  be  proven  from  the  Fourier  transform  definition. 

Equation  (3.12)  is  the  basis  for  the  rest  of  the  work  in  this  report.  (With 
some  work,  one  can  show  its  equivalence  to  equation  (2.50)  in  Hunt  and 
Graham  (1978);  this  is  demonstrated  in  app  B.)  Its  main  utility  is  that  it 
gives  the  vertically  inhomogeneous,  boundary-influenced  cross  spectra  (j)ij 
entirely  in  terms  of  the  homogeneous  cross  spectra  ^if\  which  were  given  in 
section  2.6.  For  example,  the  autospectrum  {i  -  j,  z  =  z')  for  a  component 
of  the  horizontal  velocity  is 

(j>n{Ki,K2\z,z)  =  4u^('^i>'^2;0)  — -^e~^^^4>^-^\Ki,K2\z) 

-h  f — )  e”^'''*^4f^(Ki,K2;0).  (3-14) 

\«/l/ 

(T4\ 

(In  deriving  this  result,  use  was  made  of  the  fact  that  <^13  '  is  purely  imagi¬ 
nary;  see  sect.  2.6.)  For  the  vertical  velocity  autospectrum,  we  have 
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(3.15) 


033(«1,K2;^,2)  =  (l  +  e  4?^(ki,K2;0) -2e 

Equation  (3.15)  is  the  same  as  equation  (2.12a)  in  Hunt  (1984),  after  one 
performs  the  integration  over  K2  in  Hunt’s  equation. 

The  first  term  in  equation  (3.12)  is  the  spectrum  that  would  be  observed 
if  there  were  no  boundary  effects.  The  remaining  terms  decay  exponentially 
with  increasing  height  and  decreasing  horizontal  eddy  scale. 


3.2  Modification  of  One-Dimensional  Correlations  and 
Cross  Spectra  Near  a  Boundary 


Unlike  the  homogeneous  case,  in  general  it  is  not  possible  to  find  correla¬ 
tion  functions  or  ID  cross  spectra  from  (f>ij  analytically.  One  must  compute 
the  Fourier  transform  or  integration  numerically.  For  example,  the  ID  cross 
spectra  in  the  x-direction  are  given  by 

/CX> 

4>ij{Ki,K2;z,z')dK2.  (3.16) 

-OO 


We  can  find  the  ID  correlation  function  i?ij(r,  0,0)  by  numerically  calculat¬ 
ing  the  inverse  Fourier  transform  of  ©ij(K,  0,0).  This  is  not  difficult  with  a 
fast  Fourier  transform  (FFT). 

Some  example  results  are  shown  in  figures  4  and  5.  The  calculations  are  for 
three  cases:  2:  =  z'  =  0.01^,  z  =  0.1^,  and  the  limiting  homogeneous 

case  {z  =  00).  The  most  obvious  feature  of  the  correlations  and  spectra 

is  that  the  vertical  velocity  is  severely  attenuated  near  the  boundary,  by  a 
loss  of  long- wavelength  (low- wavenumber)  energy. 

Another  interesting  feature  is  that  Rii{r^0]z,z)  ^  R22{r,0]z,z)  near  the 
boundary  (and  hence  also  0ii(«;  0, 2;,  2:)  ©22(^;  0,  2:,  ^)).  Note  that  this 

is  not  the  case  away  from  the  boundary.  Let  us  consider  the  behavior  of 
the  horizontal  velocity  spectra  at  small  heights  in  more  detail.  Making  the 
approximations  zji  1  and  k\iZ  1  in  equation  (2.38),  we  find 


This  is  the  same  as  equation  (2.42),  except  that  the  occurrence  of  K2  in 
the  second  term  in  the  square  brackets  has  been  replaced  by  Since  (j)ii 
depends  only  on  near  the  ground,  so  too  must  ^22?  by  horizontal  isotropy. 
Hence,  when  the  2D  spectra  are  integrated,  we  find  ©n  ©22* 


16 


Figure  4:  Inhomogeneous,  autocorrelation  function  with  displacement  r  along 
xi  axis,  (a)  z  =  1  velocity  component,  (b)  z  =  2.  (c)  z  =  3. 
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Normalized  wavenumber,  k€ 


Figure  5:  Inhomogeneous,  autospectral  function,  as  a  function  of  wavenum 
her  along  ki  axis,  (a)  i  =  1  velocity  component,  (b)  i  =  2.  (c)  i  =  3. 


3.3  Parameter  Selection 


Besides  u,  which  is  normally  set  to  1/3,  the  only  two  adjustable  parameters 
in  the  correlation/spectral  model  are  the  variance  <t^  and  the  length  scale 
£.  We  need  to  select  values  for  these  parameters  that  are  appropriate  for  a 
CBL.  Caughey  and  Palmer  (1979)  found  the  variance  of  all  three  velocity 
components,  away  from  the  influence  of  the  ground,  to  be 

<7‘^  =  0.Zbwl  (3.18) 

where  w*  =  {Qgzi/To)^/^  is  the  CBL  velocity  scale,  with  Q  being  the  surface 
heat  flux,  g  the  gravitational  acceleration,  Zi  the  inversion  height,  and  To  the 
surface  temperature. 

The  length  scale  £  can  be  chosen  in  a  variety  of  ways.  I  follow  Peltier  et 
al.  (1996)  by  choosing  it  to  match  the  inertial  subrange  asymptote.  In  the 
inertial  subrange,  we  must  have 

f{K)  =  (ai/2)€2/3«-V3,  (3.19) 

where  /(k)  is  the  ID  longitudinal  spectrum,  e  the  TKE  dissipation  rate,  and 
ai  is  a  constant,  approximately  equal  to  0.52  (Kaimal  et  al.  1972).  (The  co¬ 
efficient  in  equation  (3.19)  is  Q:i/2,  rather  than  oi  as  in  Kaimal  et  al.  (1972), 
because  /(k)  is  a  two-sided  spectrum,  whereas  Kaimal  et  al.  dealt  with  a 
one-sided  spectrum.  This  distinction  is  discussed  at  the  end  of  sect.  2.5.) 
Taking  the  high  wavenumber  limit  of  equation  (2.17)  and  comparing  it  with 
equation  (3.19),  we  have  the  requirement 

f-  2r(5/6) 

Lx/^r(l/3)ai 

For  freely  convective  conditions,  e  «  0.8w^/zi  (Caughey  and  Palmer,  1979). 
Hence,  equations  (3.18)  and  (3.20)  together  imply 

£  =  0.2Zzi.  (3.21) 


(3.20) 
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4.  Comparisons  With  Previous  Results 


This  section  compares  the  model  developed  above  with  field  measurements, 
large-eddy  simulation  (LES)  data,  and  other  models  of  CBL  structure.  Be¬ 
cause  is  not  easily  determined  directly  from  field  measure¬ 

ments,  the  comparisons  I  make  are  for  a  variety  of  limiting  cases,  chosen  to 
represent  a  variety  of  vertical  and  horizontal  structural  characteristics  and 
a  broad  range  of  spectral  scales. 


4.1  Energy  Spectra 

Peltier  et  al.  (1996)  recently  developed  a  model  for  horizontal  spectra  in 
the  atmospheric  boundary  layer.  Their  model,  which  included  the  surface 
blocking  effect  on  large  eddies,  agreed  quite  well  with  measured  spectra  for 
both  the  horizontal  and  vertical  velocities. 

Peltier  et  al.  developed  equations  for  the  TKE  of  the  horizontal  and  vertical 
velocities,  defined  in  my  notation  as 

EhiKh-,z)  =  TTK/i  [  4>n{Ki,  K2;  z,  z)  +  <f)22{Ki,  K2-,  z,  z)] ,  (4.1) 


and 


Ev{Kh-,  z)  =  -KKh(p23{Kl,  K2\ Z,  z). 


Their  model  for  Eh  was 


E^(Kh)  = 


Ci£iS^Kh  _  CpK'h^p 
[C2  +  (k/,4)2]4/3  (1  +  «;|^2)4/3  ’ 


(4.2) 


(4.3) 


where,  in  the  second  version,  I  have  simplified  the  equation  using  the  sub¬ 
stitutions  ip  =  ill  and  Cp  =  c\s^  jc2^ ,  The  superscript  p  indicates  the 
Peltier  et  al  model.  The  vertical  TKE  spectrum  is  related  to  the  horizontal 
spectrum  by  a  transfer  function  T{Kh^  z): 


E^{Kh)=^T{Kh,z)El{Kh)^  (4.4) 


The  transfer  function  is  tailored  “to  meet  the  continuity  constraint  [at  the 
lower  boundary]  and  to  maintain  the  proper  level  in  the  inertial  subrange.” 
Peltier  et  al  select  the  constants  and  the  transfer  function  such  that 


and 


Eh{Kh-,z) 


O.SbKhzf  2 
[23  +  (/c,Zi)f/3“'*’ 


(4.5) 


T{Kh-,z) 


{KhZ)'^/2 
0.62  +  |(k/jz)2' 


(4.6) 
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The  equations  for  Eh  and  Ey  corresponding  to  this  report’s  model  follow 
from  equations  (3.14)  and  (3.15),  after  the  appropriate  spectral  forms  are 
substituted  from  section  2.6.  The  results  are 


Eh{Kh;z) 


(1  +  k^£2).+i 


- _=e-'‘'>^X.+i(a) 

2''T(v')\  I  +  kI^ 


(4.7) 


and 


Ev{Kh\z) 


a'^Kle^ 

(1  +  Kipy+^ 


[i^iu  +  1)  (l  +  e-2"'*^) 


2'^r(z/) 


e-'''*^A',.+2(a) 


(4.8) 


The  two  models  are  compared  in  figure  6,  for  2  =  10  m  and  Zi  =  1000  m. 
The  differences  between  the  models  are  quite  small,  well  within  the  scatter  of 
the  field  measurements  plotted  in  Peltier  et  al  (1996),  Both  models  predict 
strong  damping  of  large-scale  (small  vertical  velocity  structure  near 

the  surface. 


4.2  Variances 

By  setting  n  =  r2  =  0  in  equation  (2.4),  one  has 

/oo  poo 

/  (j)ij{Ki,K2;Z,Z^)dK.idK2^  (4.9) 

-cx>  J—oo 


Figure  6:  Comparison  of  energy  spectral  models  from  Peltier  et  al.  (1996) 
and  this  work.  A  height  of  z  =  lOcalculations. 
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The  height-dependent  (co)variances  are  defined  as 

afj{z)  =  Rij{0,0,z,z).  (4.10) 

In  general,  we  must  perform  the  2D  integration  in  equation  (4.9)  numerically. 
Alternatively,  we  can  integrate  the  energy  spectra  in  one  dimension,  e.g., 

roo 

o-33(z)  =  2/  EviKfi]  z)  dKfi.  (4.11) 

Jo 

Model  predictions  for  the  variances  are  plotted  in  figure  7,  and  compared  to 
data  from  Caughey  and  Palmer  (1979)  and  an  LES  of  a  highly  convective 
CBL.  Caughey  and  Palmer’s  data  are  a  combined  set  from  the  Minnesota 
and  Ashchurch  experiments.  The  LES  data  were  supplied  by  S.  Khanna  and 
J.  Brasseur  of  the  Pennsylvania  State  University.  Their  code  was  based  on 
Moeng’s  (1984),  and  the  computational  grid  for  the  run  was  5  km  on  a  side 
in  the  horizontal  and  2  km  in  the  vertical,  with  128^  grid  points.  The  surface 
heat  flux  was  0.24  K-m/s,  and  the  geostrophic  wind  speed  was  1  m/s.  The 
resulting  value  for  —ZijL  (where  L  is  the  Monin-Obukhov  length)  was  730. 

Although  there  is  much  scatter  to  the  field  data,  there  is  generally  good 
agreement  with  the  model  predictions.  In  particular,  the  near-surface  height 
dependence  of  the  vertical  velocity  variance  is  predicted  extremely  well.  The 
variances  for  the  LES  data  tend  to  be  higher  than  either  the  predictions  or 
field  data,  by  roughly  25  percent. 

Let  us  analyze  the  near-surface  height  dependence  of  0-33(2)  in  detail.  When 
2  -C  most  of  the  spectral  contribution  to  0-33(2)  comes  from  wavenumbers 


Figure  7:  Velocity  variances  as  a  function  of  height. 


such  that  KfiZ  <C  1.  Hence,  in  equation  (3.15),  we  can  make  the  approxima¬ 
tion  exp(— Ac/iz)  «  1,  resulting  in 

<^33(^1,  K2;  2,  2)  =  K2‘,0)  —  2<f)^^\Kl,K2',z).  (4.12) 

Taking  the  Fourier  transform,  and  setting  ri  =  r2  =  0,  we  have 

a|3(2)  =  2a2-24f  (0,0,2).  (4.13) 


The  function  rI^\o,0,z)  is  an  example  of  a  longitudinal  autocorrelation 
function:  the  displacement  is  taken  in  the  same  direction  as  the  velocity 
component.  The  longitudinal  autocorrelation  /(r)  =  R^^\o,Q,r)/a‘^  was 
given  earlier  as  equation  (2.18).  Expanding  the  modified  Bessel  function  for 
small  arguments,  we  find 


0-33(2) 


2a^Til-u)  fzy' 
i/r(z/)  v2^y 


(4.14) 


Setting  the  parameters  i/,  <7^,  and  £  as  discussed  in  section  2,  we  arrive  at 

cUz)  «  l.sv,l  (4.15) 

This  result  is  identical  to  Caughey  and  Palmer’s  (1979)  equation  3.  Since 
Wfc/w*  =  (z/zi)^^^  (where  Wfc  is  the  local  free  convection  velocity  scale), 
the  result  agrees  with  previous  findings  that  the  vertical  velocity  obeys  local 
free  convective  scaling  near  the  surface  (Wyngaard,  Cote,  and  Izumi,  1971). 
In  the  model  in  this  report,  the  free  convective  scaling  has  appeared  as  a 
natural  consequence  of  the  3D  correlation/spectral  model. 


4.3  Vertical  Velocity  Correlations 


In  this  section,  I  consider  correlations  of  the  vertical  velocity  at  two  different 
heights,  1133(0,0,2,2').  This  function  can  be  evaluated  from  equation  (4.9). 
It  is  convenient  to  convert  the  correlation  function  to  a  correlation  coefficient, 
by  normalizing  by  the  variance  at  each  height: 


P33{z,z') 


i?33(0>  ^-.Z^z') 
(zzz{z)azz{z')  ' 


(4.16) 


Hunt  et  cd.  (1988)  considered  models  for  ^33(2, 2').  Based  on  a  heuristic 
argument  involving  image  eddies,  they  concluded  that 

z 

for  z-^z'  <^Zi.  Hunt  et  at  found  that  this  simple  equation  agrees  quite  well 
with  experimental  data,  as  shown  in  figure  8.  Also  shown  on  the  plot  are 
predictions  from  the  present  report’s  model  for  several  values  of  z'  jzi.  There 
is  good  agreement  when  z'  is  less  than  about  0.l2i,  but  as  2'  is  increased, 
the  models  diverge.  This  is  not  particularly  surprising,  since  Hunt  et  aVs 
equation  strictly  applies  only  for  small  z'  jzi.  Hunt  et  aVs  experimental  data 
cover  a  range  of  z/zi  values  from  about  0.01  and  0.32.  Most  of  their  data 
points  for  small  z/z'  correspond  to  relatively  large  values  for  z'/zi- 
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Figure  8:  Correlation  coefficient  for  vertical  velocity,  solid  line 

represents  Hunt  et  aVs  (1988)  equation.  Present  model  is  represented  by 
dashed  lines,  for  which  three  values  of  2'  are  shown. 

4.4  One-Dimensional  Spectra 

Longitudinal  autospectra  ©Qa(^i;  0,  z^)^  for  each  of  the  three  velocity  com- 

ponents,  are  plotted  in  figure  9  and  compared  to  data  from  the  Minnesota 
experiment  (Kaimal,  1978).  The  ID  spectra  were  determined  numerically 
by  integration  of  2D  spectra.  Predictions  are  shown  for  two  heights,  2:  =  4 
and  jz;  =  32  m,  as  Kaimal’s  measurements  included  data  taken  from  vari¬ 
ous  hdghts  between  these  limiting  values.  The  inversion  height  was  taken  as 
Zi  =  1000  m  for  the  predictions.  (When  the  predictions  are  normalized  as 
in  figure  9,  they  are  not  very  sensitive  to  the  value  of  Zi.)  The  spectra  for 
the  horizontal  velocity  components  are  normalized  by  the  variance,  which  is 
height  dependent.  The  vertical  velocity  spectrum  is  normalized  by 
where  0.75\z/L\‘^^^  for  convective  conditions  (Kaimal,  1978). 

Agreement  between  the  measured  spectra  and  the  predictions  is  fairly  good. 
For  the  horizontal  velocities,  there  is  more  high-frequency  energy  in  the 
measurements;  shear-generated  turbulence  is  the  likely  cause.  The  cross-wind 
component  of  the  horizontal  velocity  (©22)  has  slightly  more  low-frequency 
energy  than  the  predictions.  For  the  vertical  velocity,  the  spectral  peak  in 
the  prediction  tends  to  be  sharper  than  the  data. 
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ID  autospectrum  Ke33/«24.“  1 D  autospectrum  KGjj/a  ID  autospectmm 


Figure  9:  One-dimensional,  longitudinal  spectral  densities:  (a)  u-component; 
(b)  u-component;  (c)  m-component. 


4.5  Vertical  Cross  Coherence 


The  vertical  cross  coherence  for  the  ui  velocity  component  is  defined  as 


Cohn  (/Cl) 


efi(/ci,0,  z,z') 

©ii(ki,  0, 2, 2)©ii(/{i,  0, 2',  z')  ■ 


(4.18) 


Cross  coherences  such  as  equation  (4.18)  often  have  been  modeled  in  the 
literature  by  the  use  of  Davenport  similarity  (Davenport,  1961).  Davenport’s 
two  hypotheses  are  (1)  the  coherence  depends  only  on  An  =  fAz/V  = 
KiAzf27T^  where  Az  =  z^—z]  and  (2)  the  dependence  on  An  is  approximately 
exponential,  Coh(/ci)  =  exp(-’aAn),  where  a  is  a  constant  or  depends  in  a 
simple  manner  on  height,  stability,  etc  (Soucy,  Woodward,  and  Panofsky, 
1982).  In  this  section,  I  ascertain  the  extent  to  which  the  spectral  model 
discussed  in  this  report  satisfies  Davenport  similarity,  and  also  compare  it 
to  experimental  measurements  of  the  vertical  cross  coherence. 

We  compute  the  coherence  Cohn  by  integrating  the  2D  cross  spectrum,  as 
in  equation  (3.16).  Plots  are  shown  in  figure  10.  For  comparison,  I  also  plot 
experimental  data  recorded  at  White  Sands  Missile  Range,  NM  (Panofsky 
and  Dutton,  1984),  and  the  Davenport  exponential  model,  exp(— aAn)  with 
a  =  10.  The  value  a  10  is  thought  to  be  representative  of  moderately  con¬ 
vective  conditions  (Soucy,  Woodward,  and  Panofsky,  1982).  Both  the  spectral 
model  and  the  Davenport  model  agree  fairly  well  with  the  data.  There  is  a 
trend  in  the  data  for  greater  coherence  at  large  An  than  is  predicted  by  the 
spectral  model. 


Figure  10:  Vertical  coherence  of  longitudinal  velocity,  as  a  function  of  An  = 
KiAzl27r.  Model  predictions  are  shown  for  three  different  combinations  of  2 
and  2'. 


4.6  Integral  Length  Scales 


The  integral  length  scale  in  the  xi  direction  is  defined  as 

1 

=  -TTT  / 

^  <yij{z)  Jo 

and  similarly  for  Lijfi{z).  By  Fourier  transforming  equation  (2.6)  with  re¬ 
spect  to  K2,  and  then  evaluating  the  result  at  ki  =  0,  r2  =  0,  one  can  show 

TT 

(2)  =  -T7~\  /  (O’  ^2;  2,  z)  dK2-  (4-20) 

•'-00 

Let  us  take  a  moment  to  consider  the  limiting  values  for  the  length  scales  far 
from  the  ground’s  influence  {z  £)  and  close  to  the  ground  {z  i).  In  the 
former  case,  we  can  use  the  spectra  developed  for  homogeneous,  isotropic 
turbulence.  It  can  be  shown,  by  integrating  equation  (2.42),  that 

iW  ^  V^r(i/  +  1/2)^  ^  ^  ^  Q  21) 

r(i/) 

Furthermore,  since  (0, ki ; 0)  =  <p[^\ku0-,0)  by  isotropy,  integration 

yields 

^(H)  ^  ^r{u+  l/2)^  ^  ^  0.086zi.  (4.22) 

2T{u) 

The  quantity  l[^\  is  an  example  of  a  parallel  length  scale  Ly,  in  which  the 
spatial  integration  is  performed  in  the  same  direction  as  the  velocity  axis.  On 
the  other  hand,  is  a  perpendicular  length  scale  L±^  since  the  integration 

is  normal  to  the  velocity  axis.  ^3^1,  of  course,  is  also  a  perpendicular  scale. 
In  homogeneous,  isotropic  turbulence,  Ly  =  2L±  (Batchelor,  1953). 

At  the  ground  level,  on  the  other  hand,  it  can  be  shown  from  equation  (3.17) 
that 

in,i(0)  =  i22,i(0)  =  »=  O.lUi. 

{t')  (4.23) 

It  is  interesting  that  near  the  ground,  the  parallel  and  perpendicular  length 
scales  become  equal,  and  are  equal  to  two-thirds  times  the  free-space  parallel 
scale.  The  length  scale  for  the  vertical  velocity  vanishes  at  ground  level,  as 
is  evident  from  equation  (3.15). 

Model  predictions,  field  data,  and  LES  data  for  the  integral  length  scales  are 
compared  in  figure  11.  The  LES  run  was  described  in  section  4.2.  The  field 
data  were  recorded  during  the  Air  Mass  Transformation  Experiment  (AM- 
TEX),  and  analyzed  by  Lenschow  and  Stankov  (1986).  Note  that  Lenschow 
and  Stankov  found  little  azimuthal  dependence  for  the  horizontal  length 
scales,  and  hence  combined  the  data  for  the  two  horizontal  velocity  compo¬ 
nents. 
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Figure  11:  Integral  length  scales  from  AMTEX  (Lenschow  and  Stankov, 
1986)  (circles),  an  LES  of  the  CBL  (dashed  lines),  and  this  report’s  model 
(solid  line).  Left:  two  horizontal  velocity  components;  right:  vertical  velocity. 


On  the  whole,  agreement  between  the  integral  length  scale  predications  and 
the  data  is  somewhat  poorer  than  for  the  quantities  considered  in  previous 
sections.  In  particular,  the  length  scales  from  AMTEX  tend  to  be  several 
times  larger  than  the  model  predictions.  The  LES  data  generally  fall  some¬ 
where  between.  But  there  is  good  qualitative  agreement,  at  least,  as  to  the 
height  dependence  of  the  length  scales  for  the  different  velocity  components. 
The  length  scale  for  the  vertical  velocity  L334  starts  out  small  near  the 
ground  and  increases  with  height.  The  model  and  LES  both  have  1/22,1  de¬ 
creasing  with  height,  although  the  surface  value  is  about  twice  as  large  for  the 
LES.  Similarly,  both  the  model  and  LES  have  Ln^  decreasing  with  height. 
However,  in  the  LES  results,  Lu,!  nearly  vanishes  at  the  surface.  This  could 
result  from  the  finite  resolution  of  the  LES. 
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5.  Conclusion 


A  3D  correlation/spectral  model  for  turbulent  velocities  in  the  CBL  was 
developed  in  this  report  and  compared  with  data  for  a  number  of  limit¬ 
ing  cases,  such  as  variances,  vertical  correlations,  ID  spectra,  and  integral 
length  scales.  In  most  cases,  agreement  was  quite  good.  Although  previous 
researchers  have  developed  simpler  formulations  than  the  3D  model,  which 
also  compare  favorably  with  the  data  and  are  valid  for  some  of  the  individual 
cases,  the  3D  model  is  based  on  a  single,  unified  approach  to  CBL  statistics. 
And  most  importantly,  it  can  be  used  in  situations  where  complete  informa¬ 
tion  on  3D  turbulence  statistics  is  required.  Future  work  on  the  model  could 
include  incorporation  of  scalar  fields  (e.g.,  temperature  and  water  vapor) 
and  incorporation  of  ground  roughness  effects. 
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App6ndix  A.  More  on  TKE  Spectrum  and  Production 


As  discussed  in  the  introduction  to  the  main  body  of  this  report,  the  scales 
of  motion  in  a  turbulent  flow  can  be  roughly  partitioned  into  three  sub¬ 
ranges:  energy-containing,  inertial,  and  dissipative.  The  energy-containing 
range  consists  of  spatial  scales  A  (~  1/k)  much  larger  than  i  (and  hence 
At  <C  1/^))  the  inertial  subrange  of  scales  smaller  than  I  but  larger  than  rj 
{If I  K  l/v):  the  dissipative  subrange  of  scales  smaller  than  r) 

(k  '>  I Iri).  The  length  scale  77  is  called  the  Kolmogorov  microscale,  and  it  is 
equal  to  where  1/  is  the  kinematic  viscosity  and  e  the  dissipation 

rate  of  turbulent  kinetic  energy  (TKE)  (Kolmogorov,  1941). 


The  TKE  model  used  in  this  report,  equation  (2.14),  applies  only  to  the 
energy-containing  and  inertial  subranges.  This  is  not  normally  a  problem  in 
acoustics,  unless  one  is  interested  in  ultrasonic  frequencies.  But  in  order  to 
tie  up  some  loose  ends  and  explain  how  figure  1  was  created,  I  discuss  in  this 
appendix  how  equation  (2.14)  can  be  extended  to  the  dissipative  subrange. 

The  starting  point  is  the  following  equation,  due  to  Corrsin  (1964),  which  is 
valid  only  for  the  inertial  and  dissipative  subranges: 


E{k)  =  exp 


(A-1) 


where  a  is  a  constant  approximately  equal  to  1.6.  By  approximating  equation 
(A-1)  to  small  kt?,  equation  (2.14)  in  the  body  of  the  report  to  large  kL, 
and  equating  the  results,  we  have  the  following  constraint  for  the  inertial 
subrange: 


Substituting 
we  find 


/3  _  4r(17/6)  2g-2/3 


(A-2) 


this  result  into  equation  (2.14),  and  combining  with  (eq  A-1), 


kH^ 


E{k)  =  (1^^2^2)1776 


--a(K?7)'‘/® 


(A-3) 


This  equation  should  work  reasonably  well  throughout  the  turbulent  spec¬ 
trum.  It  is  plotted  in  figure  1(a)  in  the  body  of  the  report  for  i  =  IO'^tj. 

Tennekes  and  Lumley  (1972)  show  that  the  flux  of  TKE  through  the  wavenum- 
ber  K  is 


T{k)  =  a-3/2Ac5/2£;3/2.  (A-4) 


Using  equation  (A-3),  we  find 


--a{Krjfl^ 


(A-5) 


Note  that  in  the  inertial  subrange,  T{k)  -  e.  This  is  because  energy  is  neither 
produced  nor  destroyed  in  the  inertial  subrange;  it  is  merely  transferred  to 
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smaller  scales,  eventually  being  dissipated  by  viscosity.  The  production  of 
TKE  at  a  given  wavenumber  is 


P{k) 


dTjn) 

dn 


(A-6) 


TKE  production  using  equation  (A-5)  was  plotted  in  figure  1(b).  Note  that 
production  is  positive  in  the  energy-containing  subrange,  and  negative  in 
the  dissipative  subrange,  as  expected.  The  total  area  under  the  curve  sums 
to  zero,  since  the  net  production  of  TKE  must  be  zero  if  conditions  are 
stationary. 


Appendix  B.  Derivation  of  the  Two-Dimensional  Cross 

Spectral  Equation  Prom  Hunt  and  Gra¬ 
ham’s  (1978)  Result 


Hunt  and  Graham  (1978)  considered  the  spectra  for  free-stream  turbulence 
above  a  plane  boundary.  The  turbulence  was  generated  by  a  flow  moving 
through  a  grid,  with  the  boundary  moving  at  the  same  speed  as  the  mean 
flow.  Hence,  there  was  no  shear  generation  of  turbulence  at  the  boundary, 
and  the  only  effect  of  the  boundary  was  to  block  the  flow.  This  makes  Hunt 
and  Graham’s  spectral  model  appropriate  for  the  convective  boundary  layer 
(CBL),  particularly  when  the  turbulence  generated  by  the  mean  wind  is  neg¬ 
ligible  in  comparison  to  convectively  generated  turbulence.  The  assumptions 
that  go  into  the  Hunt  and  Graham  model  are  actually  equivalent  to  those 
used  in  this  report  (sect.  3);  mainly,  the  blocking  flow  at  the  boundary  is 
assumed  to  be  irrotational.  Hence,  Hunt  and  Graham’s  final  results  should 
be  equivalent  to  ours.  The  equivalence  is  demonstrated  in  this  appendix. 

Hunt  and  Graham’s  equation  (2.50),  for  the  one-dimensional  (ID),  longi¬ 
tudinal  cross  spectra,  is  the  starting  point  for  this  proof.  Rewriting  that 
equation  so  that  the  coordinate  system  follows  the  atmospheric  convention 
{z  vertical),  we  have 

/cx>  roo 

/  (x,  2;  K)Mjmix  ,z;k) 

-00  J —00 

X  exp[iK2iy  —  d«3-  (B-1) 


Here  is  the  “homogeneous”  3D  spectrum:  the  spectrum  generated  by 
the  grid  that  would  persist  without  any  blocking  from  the  boundary.  The 
tensor  Mij  is  given  by  Hunt  and  Graham  as 


Mij{x,Z\K) 


0 

0 


0 

gi/caz+iKiX 

0 


(B-2) 


where  =  Equivalently,  one  Can  write 

Mij(x,  2;  k)  =  +  mi(Ki,  /C2)5j3e“'^'‘'")  ,  (B-3) 


with 


mi{Kl,K2)  =  {iKl/Kh,iK2/Kh,—l)- 


(B-4) 


The  2D  cross  spectrum  (j)ij  follows  from  ©y  if  we  omit  the  Fourier  trans¬ 
formation  integral  involving  K2-  Doing -so,  and  setting  x  =  x',  y  =  y',  we 
find  /  \ 

<^ii(Ki,«2;^,^')=  r  M;,(0,2;K)Mj>„(0,2';/c)$J^^(K)d«3. 

(B-5) 
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Substituting  with  equation  (B-3),  we  find 


/oo 

-OO  *- 


+ 

+ 

+ 

K2)Tnj{Kl,  K2)^33\k)  dKz- 

(B-6) 

The  integrations  can  be  calculated  explicitly,  with  the  Fourier  transform 
relation 

roo 

K2;  rs)  =  J  ^(k)  exp(z/C3r3)  dK^. 

(B-7) 

The  result  of  performing  these  integrations  is  equation  (3.12).  Hence,  Hunt 
and  Graham’s  result  is  indeed  equivalent  to  the  result  in  this  report. 
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